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IS CALVERA A GAMMA-RAY PULSAR? 
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ABSTRACT 

Originally selected as a neutron star (NS) candidate in the ROSAT All-Sky Survey, 1RXS 
J141256.0+ 7 92204 ("Calvera'') was discovered to be a 59 ms X-ray pulsar in a pair of XMM-Newton 
observations (jZane et al.l 120111 ). Surprisingly, their claimed detection of this pulsar in Fermi 7-ray 
data requires no period derivative, severely restricting its dipole magnetic field strength, spin-down 
luminosity, and distance to small values. This implies that the cooling age of Calvera is much younger 
than its characteristic spin-down age. If so, it could be a mildly recycled pulsar, or the first "or- 
phaned" central compact object (CCO). Here we show that the published Fermi ephemeris fails to 
align the pulse phases of the two X-ray observations with each other, which indicates that the Fermi 
detection is almost certainly spurious. Analysis of additional Fermi data also does not confirm the 
7-ray detection. This leaves the spin-down rate of Calvera less constrained, and its place among the 
families of NSs uncertain. It could still be either a normal pulsar, a mildly recycled pulsar, or an 
orphaned CCO. 

Subject headings: pulsars: individual (1RXS J141256.0+792204, PSR J1412+7922, Calvera) — stars: 
neutron 



1. INTRODUCTION 

The NS c andidate 1RXS J1412 56.0+792204, dubbed 
"Calvera" ([Rutledge et al.l [2008h . was selected from 
the ROSAT All- Sky Survey, and observed by Chan- 
dra (|Rutledge et all l20ul (Shevchuk etUI [2001 . A 
deep radio pulsar s earch showed that it is radio quiet 
([Hessels et al.ll2007f) . It was not until a pair of XMM- 
Newton obs ervations was obta ined with high time res- 
olution that iZane et all ([201 If ) discovered 59 ms pulsa- 
tions from Calvera. The classification of Calvera among 
the families of NSs is not yet understood. Its X-ray 
spectrum is characterized as a blackbody of tempera- 
ture s s 0.2 keV, or a hydro g en atmosphere of T f» 0.1 
keV ([Shevchuk et all [2001 IZane et al.l I20T1 . Two- 
temperature models provide a better fit, and the sur- 
face temperature must be nonuniform because pulsa- 
tions are seen. Calvera's properties di stinguish it from 
seven isolated NSs (INSs: Habcrl 2007), also discovered 
by ROSAT, which are slowly rotating (P = 3 - 11 s), 
cooler NSs in the solar neighborhood. X-ray timing and 
spectroscopy and kinematic studies of the INSs indicates 
that they have strong m agnetic fields, B s » 2 x 10 13 G , 
and are rj 10 6 years old ([Kaplan fc van Kerkwiik| [2009). 
Calvera is at least twice as hot as the INSs, indicating 
an age < 10 5 yr accord ing to minimal NS cooling curves 
(|Paee et al.ll2004l2009D . Even though it is at high Galac- 
tic latitude (£, b) = (118°, +37°), if Calvera is a passively 
cooling NS it must be close to its birth place in the disk, 
implying a maximum distance of a few hundred parsecs. 
Depending on the X-ray spectral model fitted, the col- 
umn density is consistent with a range of values that 
does not furth er constrain the d istance. Based on the 
spectral fits of Za ne et all ([201 If ), the bolometric flux is 
uncertain by about a factor of 2, and the luminosity is 
L x « 1.7 x 10 3i 

^300 er § s X j where d^oo is the distance 
in units of 300 pc. Calvera remains radio quie t even after 
a dee per search for radio pulsations at 59 ms (jZane et al.l 

HSU). 



Analyzing d ata from t h e Fe rmi Large Area Tele- 
scope (LAT), IZane et all ([201 If) claimed that 59 ms 
pulsations are detected from 1RXS J141256.0+792204 
at > 100 MeV. Apart from the marginal significance 
of the detection, this result is surprising because their 
ephemeris requires no frequency derivative over the 21 
month time span analyzed. Their effective 2a upper limit 
is \f \ < 2.6 x 10~ 16 Hz s _1 , implying spin-down power 



E = -4:ir 2 Iff < 1.7 x 10 32 
r c = -f/2f > 1.0 x 10 9 yr, 



erg s , characteristic age 
and magnetic field strength 

B s = 3.2 x 10 19 ^-//-3 < 7.4 x 10 9 G. Here we assume 

a moment of inertia I = 10 45 g cm 2 , B s is the equa- 
torial surface dipole field, and frequency / = 16.89 Hz. 
(We are unable to account for their quoted upper limit 
B s < 5 x 10 10 G, which would allow / = -1.2 x 10~ 14 
Hz s" 1 .) 

These timing parameters imply that 1RXS 
J141256. 0+792204 is not just a passively cooling 
NS, but is converting a large fraction of its meager 
spin-down power into 7-rays. The quoted pulsed 7-ray 



luminosity is 



1.5 x 10 3 



^300 er S s assum- 



ing isotropic emission. All other 7-ray pulsars have 
E > 2 x 10 33 erg s" 1 (|Abdo et all l2010al ); Calvera 
would be the least energetic 7-ray pulsar by an order 
of magnitude. In this Letter, we show that the Fermi 
detection of Calvera is almost certainly false. We then 
discuss the implications for the nature of Calvera. 

2. XMM-NEWTON TIMING OF THE PULSATIONS FROM 
1RXS J141256.0+792204 

We reduced the now archival EPIC pn CCD data from 
the two XMM-Newton observations of Calvera listed 
in Table [TJ These were taken in small window mode 
with 5.7 ms time resolution, and were separated by 40 
days. The data were processed with SAS version xmm- 
sas_20090112_1802-8.0.0. We extracted photons in the 
0.15 — 2 keV band from a 20" radius aperture around 
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Table 1 

XMM-Newton EPIC pn Timing of 1RXS J141256.0+792204 



ObsID 


Date (UT) 


Date (MJD) 


Span (s) 


Exp. (s) 


Frequency (Hz) a 


Z\ 


0601180101 


2009 Aug 31 


55,074.30 


19,911 


13,941 


16.8924052(25) 


141.1 


0601180201 


2009 Oct 10 


55,114.18 


27,816 


19,477 


16.8924041(15) 


201.9 



1 sigma error in parenthesis. 
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Figure 1. The two XMM-Newton pulse profiles of 1RXS 
J141256.0+792204 in the 0.15 - 2 keV band, each folded on its 
peak frequency from Table [T] and aligned arbitrarily in phase. 
Background has been subtracted and the counts are normalized 
to 1. Two cycles are plotted. 

the source. Background was taken from an adjacent re- 
gion on the CCD. We applied the conversion to Barycen- 
tric Dynamical Time u sing the precise Chan dra posi- 
tion of the source from iShevchuk et al.l ([2009), R.A. = 
14 h 12 m 55?84, Deck = +79° 22' 03."7 (J2000.0). The 
0."6 position uncertainty is a negligible source of er- 
ror on the absolute timing. Table Q] lists the peak fre- 
quencies of the two XMM-Newton observa tions derived 
from a power sp ectrum (Rayleigh test; iStruttl \l88Ck 
iBuccheri et al.|[l983l) . We folded each light curve at the 
peak period. The pulse profiles are shown in Figure HJ 
where we have aligned them arbitrarily in phase. The 
two profiles have a consistent, quasi-sinusoidal or trian- 
gular shape, and a pulsed-fraction of « 18%. The precise 
agreement between the frequencies at the 10 -7 fractional 
level argues against any orbital motion (but see below), 
which is supported by the absence of an optical counter- 
part. 

As noted by iZane et akl (|2011l ) , the precision of the 
two frequency measurements is insufficient to join these 
widely spaced observations coherently and obtain more 
precise timing parameters. Since the measured frequen- 
cies agree within their errors, we can only derive an upper 
limit on the spin-down rate. Using the data in Table [IJ 
we calculate an upper limit on the frequency derivative 
by propagating the errors and dividing the difference of 
the frequencies by the time interval between the obser- 
vations. The 2tT upper limit is |/| < 2.0 x 10~ 12 Hz s _1 . 
The corresponding 2a limits on the spin-down proper- 
ties are E < 1.3 x 10 36 erg s _1 , r c > 1.3 x 10 5 yr, and 
B s < 6.5 x 10 11 G. These are the limits that we will 
conclude are the best curren t ly ava ilable. They are con- 
sistent with the IZane et akl (|2011| ) results of the same 
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Figure 2. The same d ata as in Figure [T] now folded on the sup- 
posed Fermi ephemeris ( Zane et al. 2011) that spans these epochs 
(see text). The bin size of the 10 bin light curve is 5.9 ms, compa- 
rable to the CCD frame time of 5.7 ms. Evidently the ephemeris 
does not phase-align the X-ray pulse profiles. 

analysis. 

Using Fermi LAT data, IZane et akl (|2011|) claimed to 
detect 7-ray pulsations from 1RXS J141256.0+792204by 
searching coherently a 21 month span from 2008 Au- 
gust to 2010 April using the Zf test. Although no 
imaging detection of a source at this position was re- 
ported, searching both / and / around the X-ray mea- 
sured values, they found a peak power Z\ = 26.2 with 
the following ephemeris: epoch to = 55094 MJD, / = 
16.892401975(2) Hz, / = -1.2(7) x 10~ 16 Hz s" 1 . Since 
this ephemeris spans the dates of the XMM-Newton ob- 
servations in 2009 August and October, it can be used 
fold the X-ray photons, and the pul s e profi les so derived 
should align in phase. IZane et akl (|2011f ) did not per- 
form this test. The result is shown in Figure 2. The 
test fails because the light curves folded according to the 
trial ephemeris are out of phase by 0.35 cycles (21 ms). 
As we discuss below, this is much larger than the maxi- 
mum ~ 3 ms uncertainty due to the 5.7 ms CCD frame 
time. Therefore, we conclude that the ephemeris repre- 
sents a noise peak in the power spectrum of (probably) 
background photons, not a true detection of pulsations. 

This test does not involve the phasing of the 7-ray 
pulse with respect to the X-ray pulse, nor is it sen- 
sitive to uncertainties on the Fermi ephemeris param- 
eters. This is because the coherent ephemeris is a 
phase ephemeris. It specifies the phase of the pulsar 
cj)(t) = f(t - t Q ) + (l/2)/(t - t ) 2 at all times during 
the span of the ephemeris, and the XMM-Newton ob- 
servations are contained well within that span. The 
quoted Fermi uncertainties, cr 7 (/) = 2 x 10~ 9 Hz and 



cr 7 (/) = 7 x 10~ 17 Hz s~\ 
drifts in relative phase, A(j> - 



can contribute only tiny 
: a 7 (f)T x = 0.007 cycles 
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and Acj) = 0.5ct 7 (/)T| = 0.0004 cycles, respectively, 
over the Tx = 40 day interval between the X-ray obser- 
vations. Even these are upper limits, as the errors are 
covariant. 

Our analysis does depend on the stability of EPIC 
pn timing in small window mode. Considerable effort 
has gone into calibrating the relative and absolute tim- 
ing of this particular mode and maintaining the accu- 
racy of the photon time assignments in the processing 
chain. The absolute a ccuracy of the XMM -Newton clock 
is better than 0.6 ms ([Kirsch et al.ll2004T ). We have car- 
ried out extensive investigations of pulsars using the pn 
small window mode, and empirical checks for consistency 
show that the absolute times are at least as accurate as 
the ~ 3 ms uncertainty due to the 5.7 ms CCD frame 
time. We summarize two of these studies here: A co- 
herent ephemeris for the 237 ms pulsar Geminga using 
10 XMM-Newton observations over a span of 7 years has 
phase residuals of 2 ms. We have used it to demonstrate 
the stable phase relationship between the X-ray and 
7-ray pulse of G eminga with AGILE (jPellizzoni et al.l 
120091 ) and Fermi (|Abdo et al.ll2010bl ) to < 1 ms. These 
agree with earlier results c omparing EGRET and ASCA 
(| Jackson fc Halpernll200l . In an ex tensive campaign on 
the 105 ms pulsar PSR J1852+0040 (IHafpern fc GotthelJ 
120101) . 16 XMM-Newton observations and seven Chan- 
dra observations (with time resolution 3 ms) spanning 
4.8 years are fitted by a quadratic phase ephemeris with 
rms phase residuals of 3.4 ms. This demonstrates that 
XMM-Newton and Chandra agree in absolute time to 
better than 3 ms. 

Our analysis does not make use of the EPIC MOS 
detector timi ng data that were obtained simultaneously 
with the pn. IZane et al.l (|2011l ) noted that the MOS de- 
tector's timing mode is not well calibrated, and they did 
not use it for timing analysis. 

3. DISCUSSION OF PREVIOUS FERMI ANALYSIS 

We address here reasons why one might consider the 
Fermi detection of Calvera to b e real despite our nega- 
tive evidence. First, IZane et all (|2011| ) suggest that ac- 
tually a more conservative upper limit on the frequency 
derivative from Fermi should be allowed, |/| < 10~ 15 
Hz s -1 . Still, its effect over the 40 day interval between 
X-ray observations would be negligible, and it would not 
change the outcome of our X-ray phase comparison. It is 
not clear why they entertain this possibility, since such a 
value would contribute 1.5 extra cycles of rotation over 
the 21 month span of their ephemeris. If / actually turns 
out to be —1 x 10 -15 Hz s _1 , it would be a different 
ephemeris from the published one, and the claimed detec- 
tion would be spurious. For that matter, the frequency 
of the pulsar could also turn out to differ by more than 
one Fourier bin (1/TL, where T 7 is the 21 month span 
of the Fermi data) from the published Fermi ephemeris, 
and still be consistent with the X-ray measured value. 
In this case as well, the claimed 7-ray detection is just 
noise. 

Second, it may be argued that, while the fitted 
ephemeris corresponds to the average values of / and 
/ over the time span, there could be timing noise that 
smears the pulse, while the signal is not strong enough to 
fit such trends with higher order terms. Under this hy- 



pothesis, the phase drift between the two X-ray observa- 
tions is a manifestation of timing noise. We consider this 
unlikely because the fitted / is already consistent with 
zero. Any detectable timing noise would vary the sign of 
/, which has not been seen in any isolated pulsar apart 
from glitch discontinuities. A phase drift of 0.35 cycles 
over 40 days would require an effective / = —5.9 x 10~ 14 
Hz s _1 over this time, almost 500 times the mean value 
of —1.2 x 10~ 16 Hz s _1 . This seems unlikely, as does 
a glitch that is not also detected in the 7-ray timing. 
No pulsar with T r > 2 Myr has been observed to glitch 
(jEspinoza et al.ll2011[ ). 

As an alternative to timing noise, it may be hypothe- 
sized that the X-ray phase shift over 40 days is evidence 
of binary motion, i.e, a planetary companion. Such an 
explanation would require the orbital period to be much 
less than the 21 month span of the Fermi ephemeris, 
but longer than the durations of the individual XMM- 
Newton pointings, which are 5.5 hours and 7.7 hours, 
respectively. We consider that 1 day < P or b < 100 days 
covers the applicable range. Under the binary hypothe- 
sis, the projected radius a ns sini of the NS orbit falls in 
a narrow range. The Roemer delay 2 a ns sinz/c must be 
> 21 ms to produce the phase shift, but < 59 ms so as not 
to smear out the supposed 7-ray pulsations. Combining 
these requirements, for m ns — 1.4 Mq we find that 

This corresponds to a minimum planet mass of m p = 
20/sinz Earth masses, and a maximum of 4/sinj Jupiter 
masses, the latter for P or b = 1 day. More likely the 7-ray 
ephemeris is spurious, and the 0.35 cycle phase shift is a 
random number, not evi dence of a p l anet. 

Other aspects of the IZane et aLl (|2011f ) analysis are 
unusual. By their own description, the 7-ray signal is 
marginal, with Z\ = 26.2, and its significance depends on 
the number of independent trials in the search. The trials 
can be assessed using the X-ray uncertainties on the tim- 
ing parameters in Section 2, crx{f) = 1-5 x 10~ 6 Hz and 
ax(f) = 1 X 10~ 12 Hz s _1 , and the T 7 = 21 month span 
of the Fermi data. The number of independent trials 
in the two-dimensional Fermi search should be at least 
4o- x (/)T 7 « 330 for frequency, and <7 X (f)T% « 3000 
for frequency derivative. These represent a search of in- 
dependent frequencies in the ±2cr interval around the 
XMM-Newton measured /, and independent frequency 
derivatives ranging from zero to the — 2er limit, i.e., only 
negative /. If so, the expected number of noise peaks of 
power Z\ > 26.2, obtained by multiplying the 1 x 10 6 

trials by the single-trial probability, e ~ 26,2 / 2 — 2x 10~ 6 , 
is of order unity. The oversampling by a factor of 10 
that was performed further reduces the statistical signif- 
icance. The crux of their argument must be that, since 
the value of / in the discovered signal is consistent with 
zero, almost no trials in / were needed to find it. Only 
this would allow that the chance probability of the re- 
sult is ~ 7 x 10~ 4 . But they do not display the power 
spectrum for the complete search. Instead, they say that 
they did not find false detections over a range of param- 
eters much wider than the X-ray uncertainties. It is not 
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Figure 3. Power spectrum of Fermi data covering the time span 
analyzed bv lZane et all (120111) , 2008 August 4 - 2010 April 27. The 
peak power Z\ = 24.4 is similar to their value of 26.2. The two- 
dimensional search in / and / was projected onto one dimension 
in this plot. 

stated what they consider a false detection, and we are 
not shown the values of the highest peaks in the extended 
search. Therefore, we don't know what to make of this 
argument. Finally, the absence of the 7-ray source in a 
spatial image is also worrisome, especially since the posi- 
tion is at high Galactic latitude with minimal confusing 
diffuse background. So there is no supporting evidence 
of a source at this position. 

4. ANALYSIS OF NEW FERMI DATA 

For completeness, we extracted and reduced Fermi 
data using the same e vent filtering and methods de- 
scribed by iZane et all (|2011l ). We first extracted the 
identical 631 day time span, 2008 August 4 - 2010 April 
27, recovering 2750 photons, similar to their 2518 pho- 
tons. A Zf search covering their ±3cr range of Fermi 
frequency derivative recovers the candidate peak at / = 
16.892401976, as shown in Figure [3] The peak power, 
Z\ = 24.4, is consistent with their value of 26.2. 

We then applied the same method to the full data set 
now available, comprising 4764 photons collected up to 
2011 May 19 (33 months). This increases the number 
of photons by 73%. A search of the same ephemeris 
parameters does not yield increased power at the claimed 
frequency. Rather, the peak previously seen is reduced 
to Z\ — 16.0 (Figured]), indicating that it was just noise. 
We regard this result as direct support of our inference 
that the supposed Fermi detection was not real. A wider 
search of thousands of trials in / is not meaningful, for 
the reasons discussed above. 

5. CONCLUSIONS 

In summary, the incorrect phasing of the X-ray obser- 
vations of Calvera using the claimed Fermi ephemeris led 
us to conclude that the 7-ray detection is probably spu- 
rious. Then, extending the Fermi analysis from 21 to 33 
months rendered the candidate signal insignificant. Here 
we comment on the implications for the nature of Calvera 
of the less restrictive 2cr limit |/| < 2.0 x 10~ 12 Hz s~ 4 
from X-ray timing. The corresponding limits on spin- 
down properties, E < 1.3xl0 36 ergs -1 , r c > 1.3xl0 5 yr, 



Figure 4. Power spectrum of Fermi data covering an extended 
time span, 2008 August 4 — 2011 May 19, showing that the peak 
in Figure [3] is no longer significant. 

and B s < 6.5 x 10 11 G, allow three scenarios. First, 
Calvera could be an ordinary pulsar with magnetic field 
strength 1 .5a or more below the mean of the birth dis- 
tribution (jFaucher-Giguere fe Kaspil 120061 ) . Second, it 
could be a mildly recycled pulsar, formerly in a binary 
system with a high-mass companion that has since un- 
dergone a supernova explosion. In this case, its mag- 
netic field strength would be intermediate between those 
of fully recycled (millisecond) pulsars, and ordinary pul- 
sars. Third, it could be an "orphaned CCO" as we de- 
scribe next. The second and third scenarios were also 
addressed bv lZane et al.l (|2011[ ). 

The class of central compact objects (CCOs) in 
supernova remnants comprises w 10 NSs, three of 
whi ch are detected pu l sars with P = 0.1 — 0.4 s. 
SeelHalnern fc GotthelJ IpoToL l2?)TTh and IGotthelf"et~aT1 
(120101 ) for observations and overview of related the- 
ory. The CCO pulsars have weak dipole fields, in the 
range 10 10 — 10 11 G, and negligible spin-down power in 
comparison with their bolometric X-ray luminosities of 
10 33 - 10 34 erg s" 1 . These 10 3 - 10 4 yr old NSs must rep- 
resent a significant fraction of NS births. After their host 
supernova remnants dissipate, orphaned CCOs will re- 
main in the region of (P, P) space where they were born, 
which is al so where the supposed m ildly recycled pulsars 
are found (jBelczvnski et al.l 120100 . An orphaned CCO 
would be distinguished from a single, recycled pulsar by 
its residual thermal X-ray luminosity. An orphaned CCO 
could be recognized as a thermal X-ray source, depending 
on its distance, while it is up to 10 5 — 10 6 yr old (not the 
characteristic age, which is orders of magnitude larger 
than the real age of a CCO). The known CCOs have X- 
ray temperatures in the range 0.2 — 0.4 keV. Calvera, be- 
ing cooler than this and less luminous than a young CCO 
by a order of magnitude or more, could be an evolved 
stage of a passively cooling CCO. It is possible, there- 
fore, that 1RXS J141256.0+792204 is the first orphaned 
CCO to be recognized. 

These scenarios can be distinguished by the spin-down 
rate of the pulsar. If B s — 6 x 10 11 G and r c = 
1.5 x 10 5 yr, Calvera is probably an ordinary pulsar. If 
B s = 1 x 10 11 G and r c = 4 x 10 6 yr, it could be an or- 
phaned CCO younger than r c , or a spin-powered, mildly 
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recycled pulsar. If B s = 1 x 10 10 G and E = 3 x 10 32 
erg s" 1 , its spin-down power is probably insufficient to 
heat its surface, and an orphaned CCO would be re- 
quired instead of a mildly recycled pulsar to explain its 
X-ray luminosity and temperature. Even though Calvera 
is not (yet) detected in 7-rays, an X-ray timing study 
is straightforward, and should detect its spin-down in 
w 1 yr even if its magnetic field strength is only w 10 10 G. 

We thank Eric Gotthclf for discussions and assistance 
with the data. This investigation is based on observa- 
tions obtained with XMM-Newton, an ESA science mis- 
sion with instruments and contributions directly funded 
by ESA Member States and NASA. 
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